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Abstract 

We extend the equation of state of GCG (Generalized Chaplygin Gas) to w < —1 regime and 
show, from the point of view of dynamics, that the parameters of GCG should be in the range of 
< a < 1. Also, dynamical analysis indicate that the phase w g = — 1 is a dynamical attractor and 
the equation of state of GCG approaches it from either w g > — 1 or w g < — 1 depending on the 
choice of its initial cosmic density parameter and the ratio of pressure to critical energy density. 
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1. Introduction 



CMB anisotropy 



HQS, 



Super novae 



and SDSS strongly indicate that our 



universe is spatially flat, with two thirds of the energy contents resulted from dark energy, 
a substance with negative pressure and can make the universe expanding in an accelerating 



ras 

a 



lion. Candidates for dark energy have been proposed as vacuum energy, quintessence 
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1% 14^ phantom 15] and GCGjlGj which is stemmed from the CUaplygiau 



gasHK. Present observation data constrain the the range of the equation of state of dark 







which indicates the possibility of dark energy with 



energy as —1.38 < w < —0.82 

w < — 1, debuted as Phantom 0]. The realization of w < — 1 could not be achieved by scalar 
field with positive kinetic energy and thus the negative kinetic energy is introduced although 
it violates some well known energy conditions Another important consequence of 



Phantom is the Big rip 20] or Big smash |2l( phase, in which the scale factor of the Universe 
goes to infinity at a finite cosmological time. The cosmological implications of Phantom 
have been widely studied Q, 24,3E,2i,3QQ,HE,Q,Q0Q3Q,Q 
and the Phantom model with Born-Infeld type Lagrangian has been proposed [40] and its 
generalization to brane world has been done in Ref. 41 1. 

GCG has a very simple equation of state, p„ = — m4( " +1) ; which yields an analytically 

' Pg 

solvable cosmological dynamics if Universe is GCG dominated. Another ambition of intro- 
ducing GCG is to unify unify dark energy and dark matter into one equation of state, also 
known as quartessence 42]. However, detailed numerical analysis turns out to disfavor the 
dark matter modelled by the GCG equation of state j3]. But no observation so far rule out 
the possibility of GCG as dark energy. It is quite possible for our Universe to contain a 
dark energy component modelled by GCG as well as another baryotropic fluid component 
mimicked by the equation of state p 7 = 7p 7 . Although the cosmological GCG system is 
analytically solvable when GCG is dominant, it is no longer possible when another compo- 
nent is present. We therefore resort to the phase analysis with which one can gain many 
important information without solving the dynamical equations. 

Previous study of GCG are focus on the case that w g > — 1. In fact, if we accept the 
notion that dark energy is modelled by a perfect fluid, then, current observations do not 
exclude the possibility of w g < — 1, instead, they even favor it|3]- Thus, in the parameter 
space of the GCG model, one should not exclude the regime w g < — 1 although this range of 
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equation of state can not be smoothly continued from the w g > — 1 regime. In this letter, we 
generalize the idea of GCG by considering the w g < — 1 regime of its equation of state. We 
show that the system could reach the late time de Sitter attractor from either w g > — 1 or 
w g < — 1 depending on the choice of the initial energy density and pressure. When w g < — 1, 
it will behave as phantom with a late time de Sitter attractor, and therefore won't lead to 
the catastrophic big rip. 

The equation of state of GCG has two free parameters, a and M, which could be fixed, in 
principled fitting t he mod e, to th e S„e or CMB da ta Q . However, when dealing 
with these fitting, we need to narrow down the possibilities of the range of the parameters 
to facilitate the numerical analysis. This is just part of the purpose of this current work, in 
which we constrain the parameter a from the cosmological dynamics of the system as well 
as the requirement of sound speed. Our result is in agreement with the numerical results 
obtained by other authors. 



2. Autonomous system 

A general study the phase space system of phantom scalar field in FRW universe has been 
given in Ref.j3]- For the GCG cosmological dynamical system, the corresponding equations 
of motion and Einstein equations could be written as, 



K 2 

H = - — (p^+Vi + Pg+Vg) 
p' 7 = -3F(p 7 + p 7 ) 

Pg = ~ 3H (Pg+Pg) 

H 2 = ^-(p, + p g ) (0.1) 

where k 2 = 8ttG, p 7 is the density of fluid with a baryotropic equation of state p 7 = (7 — l)p 7 , 
where < 7 < 2 is a constant that relates to the equation of state by w 7 — 7 — 1; p g is the 
energy density of GCG with an equation of state 

rg 

The over dot represents derivative with respect to cosmic time t, and H is Hubble parameter. 
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To analyze the dynamical system, we rewrite the equations with the following dimension- 
less variables 



K 2 Pg 

y - 



3H 2 

N = In a (0.3) 



The dynamical system will then reduce to 

= _3( x + y ) + 3x[7(l -x) +x + y] 

^ = 3a(y + y 2 /x) + 3y[ 7 (l - x) + x + y] (0.4) 
Accordingly, the Friedman equation yields 

tt g + fi 7 = 1 (0.5) 

2 

where Q g = x and fi 7 = |^ are the cosmic density parameters for GCG and baryotropic 
fluid respectively. The equation of state for the scalar fields could be expressed in terms of 
the new variables as 



w g = ^ = V - (0.6) 

Pg ^ 



and the sound speed is 



c 2 s = -a*- (0.7) 

x 

The critical points of the system are (x,y) = (1,0) and (1,-1) which correspond to a 
matter dominated phase and vacuum energy dominated phase respectively. If we linearize 
the system near its critical points and then translate the system to origin, we could readily 
write the first order perturbation equation as 



U' = A ■ U (0.8) 
where U is a 2-column vector consist of the perturbations of x and y. A is a 2 x 2 matrix 
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, -3 + 37 — Q'yx + 6x + 3y —3 + 3x 

A = ' (0.9) 

- ? ^-3 iy + 3y 3a + ^ + 37 -3 7 x + 3x + 6?/ 

The stability of the critical points is determined by the eigenvalues of the matrix A at the 
critical points. For the point (1,0), the two eigenvalues are 



\ x = 3-37 (0.10) 
A2 = 3a 

So, it may be stable if a < and 7 > 1. However, since we want the GCG behaves as 
dark energy, it won't be appropriate if (1, 0) corresponds to a stable attractor phase. In 
other words, a < should not be considered in the real models. While for the critical point 
(1,-1), the corresponding eigenvalues of matrix A are 



\ x = -3(1 + a) (0.11) 
A 2 = -3 7 

It is clear that (1,-1) is stable for a > — 1. This critical point corresponds to a phase that 
GCG is dominant (Q g = 1) and its equation of state is w — — 1. So, it is a late time de Sitter 
attractor |44j. Combine the above two constraints as well as the requirements of c 2 s < 1 at 
the de Sitter attractor, we can readily reach the constraint for the a parameter should be 
< a < 1. 

In the following, we study the above dynamical system numerically. For definite, we 
choose the parameters as 7 = 1 and a = 0.5. The initial x and y are chosen as shown 
in Table I and the results are contained in the following Fig.l-Fig.4. From Fig.2, one can 
observe that for different initial p g and p g , the equation of state w g could approach to the 
de Sitter phase (w g = — 1) from either w g > — 1 or w g < — 1. The critical case will be that 
the initial x and y are such chosen that w g = y/x = —1. In other words, the initial choice 
of w g > — 1 or < — 1 determines whether the equation of state of GCG will mimic that of 
quintessence or phantom. It is worth noting that if we choose the parameters so that the 
GCG behaves as phantom, it is no longer possible to make it behave as matter at early epoch 

5 



x 0.14 0.15 0.16 0.17 0.18 0.19 0.20 
y -0.19 -0.18 -0.17 -0.16 -0.15 -0.14 -0.13 

TABLE I: The initial values of x and y in the plots Fig.l-Fig.4 
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x 

FIG. 1: The phase diagram of the GCG system in terms of x and y for different initial x and y. 

because the equation of state could not smoothly evolves from to w g < — 1. However, in 
our setup of this paper, we have included a baryotropic fluid that could be used to mimic 
the matter sector of our universe and thus GCG can be considered only as dark energy. 
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FIG. 2: The evolution of the equation of state of GCG for different initial x and y. The curves from 
bottom to top correspond to the initial conditions specified in Table I from left to right respectively. 




— GCG 

— Matter 



N 

FIG. 3: The evolution of the cosmic density parameter for matter fi 7 and GCG Q, g respectively at 
different initial x and y. The plot indicates that the evolution of Q g and f2 7 is not very sensitive 
to the initial condition due to the attractor property. 
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FIG. 4: The evolution of sound speed of GCG (? s for different initial x and y. The curves from top 
to bottom correspond to the initial conditions specified in Table I from left to right respectively. 

3. Conclusion and Discussion 

In this letter, we analyze the dynamical evolution of GCG for different parameters and 
initial conditions. We show that different initial x and y will lead to different tracks (w g > 
— 1 and w g < —1) for the equation of state w g to approach the de Sitter attractor phase 
(w = —1). That is to say, the GCG could mimic both quintessence and phantom during its 
evolution depending on the initial conditions. We also give constraint to the parameter of 
the model as < a < 1 from the requirement of its dynamics and sound speed. 

On the other hand, the existing studies of GCG and its fitting to Supernovae data focus 
oiiWj > —1 2, 3]- But in fact, observations do not exclude the possibility of w g < — 1. 
So the GCG with w g < —1 should also be considered and its fitting to SNela data will be 
carried out in a preparing work. 
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